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CHAPTER  I 
INTRODUCTION 


A.  BACKGROUND 

The  analysis  of  doubly  periodic  phased  antenna  arrays  has  been  well 
documented  over  the  years.  However,  these  studies  have  generally  been 
limited  to  elements  which  lie  completely  within  the  plane  of  the  array, 
referred  to  here  as  a  totally  planar  array.  There  is  a  good  reason  for 
this.  The  nonplanar  array,  or  a  two-dimensional  array  with  elements  not 
wtihin  the  plane  of  the  array,  is  much  more  difficult  to  analyze  due  to 
the  complexity  of  the  self /mutual  impedance  terms.  The  need  for  greater 
flexibility  in  design  has  resulted  in  the  recent  study  of  arbitrarily 
oriented  dipoles  by  English  [13,  nonplanar  loop  antennas  by  Kent  [23, 

"V"  dipoles  by  Lin  [33,  and  the  Clavin  element  by  Ng  [43. 

The  difficulties  involved  can  be  seen  by  considering  a  single 
nonplanar  array  in  free  space  shown  in  Figure  1.1.  A  side  view  is  shown 
in  Figure  1.2.  English  identifies  three  regions  of  space:  Region  1 
contains  only  left  going  plane  waves.  Region  2  contains  both  left  and 
right  going  waves,  and  Region  3  contains  only  right  going  plane  waves. 
The  induced  voltage  measured  at  the  terminals  of  an  antenna  element  (2) 
due  to  the  field  generated  by  the  antenna  array  (1)  is  given  by 
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Figure  1.1.  Standard  form  of  rectangular  lattice  of  antenna  elements 
REGION  I  REGION  2  REGIONS 


Figure  1.2.  The  field  generated  by  this  nonplanar  array  takes  a 
different  form  in  each  of  these  three  regions. 


(1.1) 


^2,1 


1 


/  (R')*E^^^R')dR' 

element  2 


where  l(2)t(R')  is  the  surface  current  distribution  on  antenna  (2)  under 
transmitting  conditions,  l(2)t  is  the  current  magnitude  at  the 
terminals,  and  e(1)(R)  is  the  field  generated  by  the  array.  The  mutal 
impedance  between  the  antenna  array  (1)  and  test  element  (2)  is 

vM 

'>•2' 

When  the  test  element  (2)  and  the  reference  element  (1)  of  the  array  are 
the  same.  Equation  (1.2)  gives  the  self  impedance  of  the  phased  array. 

In  order  to  find  this  self  impedance,  the  field  next  to  the  elements,  in 
Region  2,  must  be  known.  Region  2  fields  contain  left  and  right-going 
waves  which  depend  on  both  the  source  point  and  observation  point.  The 
impedance  calculation  involves  one  integration  over  the  current  sources 
to  produce  the  field  e(1)(R)  and  another  integration  over  the  test 
element  current  distribution  l(2)t(R)  of  Equation  (1.1).  The  resulting 
expressions  for  a  nonplanar  array  involve  tedious  iterated  integrals. 
With  a  totally  planar  array,  the  integrals  are  independent  and  may  be 
separated. 

Besides  the  complexity,  these  self  impedance  equations  no  longer 
have  the  rapidly  decaying  exponential  term  which  guarantees  fast 
convergence  of  the  doubly  infinite  summation.  This  means  longer 
computer  run  time. 
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With  dielectric  layers  added,  English  showed  that  the  self 
impedance  of  an  array  may  be  decomposed  into  five  modes, 

^self  “  ^1  ^2  ^3  ^4  ^5  * 

Modes  two  thru  five  are  the  four  bounce  modes  which  exist  because  of  the 
dielectric  boundaries.  These  modes  use  the  same  pattern  factors  and 
have  the  same  form  for  nonplanar  as  for  planar  arrays.  Mode  one  is  the 
direct  mode  and  is  the  self  impedance  of  the  array  embedded  in  a 
homogeneous  space  with  dielectric  constants  cp  and  up.  This  report 
presents  a  method  for  finding  the  self  impedance  of  a  nonplanar  array  in 
a  homogeneous  space  which  can  be  combined  with  the  four  bounce  mode 
impedances  to  find  the  total  impedance  of  a  nonplanar  array  in  a 
stratified  dielectric  medium.  Dielectric  layers  will  no  longer  be 
considered  here;  all  arrays  will  reside  in  a  homogeneous  medium. 

B.  OVERVIEW 

The  typical  geometry  for  a  single  nonplanar  array,  as  shown  in 
Figure  1.1,  consists  of  a  rectangular  lattice  in  the  x-z  plane  with 
interelement  spacings  and  D^.  Skewed  grids  in  the  x-z  plane  have 
also  been  examined  by  Kornbau  [5].  As  an  example,  consider  the  array  of 

A 

y  oriented  dipoles  in  Figure  1.3  and  denote  the  self  impedance  of  this 
array  by  ZoCDj^.D^).  We  wish  to  find  this  self  impedance.  This  example 
will  show  the  intuitive  concepts  of  the  method  to  be  presented  in  the 
following  chapters  and  provide  a  first  look  at  the  specialized  notation 
invol ved. 
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Figure  1.3.  Side  view  of  y  oriented  nonplanar  dipole  array  with 
interelement  spacing  02* 

E^y  E_y 


BE  PLACED  IN  THIS  REGION 

Figure  1.4.  Volume  array  formed  by  adding  an  infinite  number  of 
identical  dipOle  arrays  with  interplanar  spacing  Dy. 
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We  begin  by  adding  an  infinite  number  of  identical  nonplanar  arrays 
with  interplanar  spacing  Dy  so  as  to  form  a  volume  array  as  shown  in 
Figure  1.4.  The  spacing  Dy  should  be  larger  than  the  total  dipole 
length.  The  field  around  our  original  two-dimensional  dipole  array  can 
be  divided  into  three  parts:  Eq  due  to  the  original  array,  E+y  due  to 
the  left  half-infinite  volume  of  dipoles,  and  E.y  due  to  the  right 
half-infinite  volume  of  dipoles.  E+y  and  E.y  are  right  going  and  left 
going  plane  waves  respectively.  This  field  may  be  integrated  over  a 
test  element  in  Equation  (1.1)  to  find  the  volume  impedance, 

^volume  “  ^  ^+y^^x’^z l*^y^  *  ^-y^^x’^z^^y^* 

(1.4) 

^o^^x’^z^  is  still  unknown.  However,  and  I  may  be  easily 

calculated  because  the  test  element  is  not  within  the  regions  which 

generate  E  and  E  .  There  is  another  way  of  calculating  Z  , 

+y  -y  j  rt  volume 

This  volume  array  is  equivalent  to  an  infinite  number  of  planar  arrays, 
all  parallel  to  the  x-y  plane,  extending  into  the  z  direction.  In  the 

A  A  A 

X-y  plane,  our  y  directed  dipoles  form  a  totally  planar  array  with  a 
self  impedance  2Q(Dj(»Dy)  which  may  be  readily  evaluated  by  standard 
methods.  Also,  Z^^  and  Z_^  can  be  calculated  to  form 

^volume  ”  ^o^^x’^y^  *  ^+z^^x’'^y l^z^  *  ^-z^^x’^y ^^z^ * 

(1.5) 
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Subtracting  Equation  (1.4)  from  Equation  (1.5)  gives 


ZJDx.D^)  -  Z„(Dx,0^)  •  Z,(Dx|D,*D^)  (1.6) 

with 

Z,(0x|D,*D^)  =  Z,^(Dx,D^|0,)  *  Z.x(Dx.Dj,|0^) 

Equation  (1.6)  shows  how  to  find  the  nonplanar  self  impedance  2^(0^, D^); 
we  calculate  the  planar  self  impedance  2^(0^, D^)  and  the  difference 
impedance  Z^(Dj(|D2'*’Dy)  and  add.  The  critical  point  in  this  procedure  is 

A 

that  the  dipole  has  maintained  the  same  y  orientation  throughout.  This 
means  that  the  vector  pattern  functions 

1  .  " 

p(r)  =  /  I(R')  e  dR'  (1.8) 

I(R  ')  element 

which  find  their  way  into  the  impedance  equations  and  which  depend  only 
upon  element  orientation,  are  the  same  in  every  impedance  term.  The 

A 

changing  lattices  have  their  effect  on  the  argument  r  and  not  on  the 
pattern  function. 

A 

For  this  example  of  a  y  directed  dipole,  we  transformed  the  planar 

A  A  A  A 

lattice  from  the  x-z  plane  in  which  the  element  is  nonplanar  to  the  x-y 
plane  where  the  element  is  planar  and  related  the  two  self  impedances 
by  Equation  (1.6).  For  an  arbitrarily  oriented  dipole  with  direction  p. 
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we  would  like  to  be  able  to  transform  our  lattice  from  the  x-z  plane  to 

A  A 

the  z-p  plane  since  this  plane  allows  us  to  easily  calculate  the  self 

A 

Impedance  of  the  totally  planar  p  directed  dipoles.  In  order  to  relate 
the  self  Impedance  of  the  dipoles  In  the  x-z  lattice  to  the  self 
impedance  In  the  z-p  lattice,  we  form  a  three-dimensional  x-z-p  volume 

A  A  A 

lattice  (x,  z,  and  p  are  assumed  to  be  linearly  independent).  This 
intermediate  volume  array  enables  us  to  write 

*  V^zl'x^P)  (1-9) 

in  which  the  terms  on  the  right  are  easily  evaluated  by  methods  which 
will  be  developed  in  Chapter  II.  Chapter  III  will  examine  this  example 
in  more  detail.  By  changing  the  lattice,  we  shall  see  how  to  find  the 
self  and  mutual  impedances  of  much  more  complex  structures. 

C.  THE  END  OF  COORDINATES 

The  plane  wave  expansion  of  the  fields  from  a  periodic  surface 
began  with  a  rectangular  lattice.  A  rectangular  coordinate  system  of 

AAA 

unit  vectors  x,y,z  was  therefore  the  natural  basis  for  expanding 
vectors.  With  the  introduction  of  skewed  lattices,  the  same  rectangular 
coordinates  were  used  with  little  difficulty  even  though  they  are  not 
the  most  natural  basis  vectors.  Using  the  method  in  this  report, 
however,  an  antenna  element  may  only  be  planar  in  a  skewed  plane  which 
does  not  contain  any  of  the  three  coordinate  axes  x,y,  or  z.  In  order 
to  calculate  the  planar  impedance  required  by  Equation  (1.6),  we  must 
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either  rotate  the  elements  into  the  x-z  plane  where  we  are  most 
comfortable  in  calculating  impedances,  or  adopt  a  different  set  of  basis 
vectors  for  the  plane  wave  expansion.  The  first  alternative  means 
giving  up  the  invariance  of  the  pattern  function  under  a  change  in 
lattice  while  the  second  means  adopting  a  more  general  theory  of 
periodic  structures.  The  choice  is  clear.  By  dropping  the  rectangular 
coordinate  system  we  are  able  to  see  directly  how  a  change  in  the 
lattice  structure  affects  the  physics  of  the  plane  wave  expansion. 
Coordinate  systems  will  only  be  used  when  it  is  necessary  to  consider 
specific  examples  in  Chapter  III, 

Chapter  II  examines  the  coordinate  free  plane  wave  expansions  in 
both  two  and  three-dimensional  lattices.  The  notation  to  be  used  for 
the  direct  and  reciprocal  lattices  is  chosen  to  remain  consistent  with 
the  classic  work  by  Leon  Brillouin  [6],  Wave  Propagation  in  Periodic 
Structures,  This  book  is  recommended  for  its  practical  insight  into  the 
mathematics  of  periodic  structures. 
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CHAPTER  II 

THE  COORDINATE  FREE  PLANE  WAVE  EXPANSION 

A.  TWO-DIMENSIONAL  ARRAYS 

The  Region  1  and  Region  3  fields  from  a  periodic  two-dimensional 
array  of  current  sources  is  derived  in  Appendix  A.  The  array  lattice, 
often  referred  to  as  the  direct  lattice,  is  defined  by  the  two  vectors 
dj  and  d2  shown  in  Figure  2.1.  The  reference  point  or  terminals  for  the 
q,m  current  source  is  located  at 

V  “  ‘’^1  "’^2  (2.1) 

We  assume  that  the  current'  at  each  lattice  point  are  identical  except 

for  a  linear  phase  difference  of  the  form  e 
gives  the  electric  field, 

Z  i(1)(r(1)) 

e“'(r)=— — —  I  i 

2iidjxd2i  n=_«B  ks-® 

(p(^^(r^)  X  r+)  X  r^ 


qiM 


.  Equation  (A. 40) 


-j6(R-R^^b«r^ 


(2.2) 
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Figure  2.1.  An  array  of  antenna  elements  defined  by  the  direct  lattice 
vectors  d^  and  62. 


with 


(2.3) 


p(l)(?)  = 


/ 


(2.4) 


element  1 


>  -(1) 

R'  is  a  vector  from  R  to  the  integration  point.  The  propagation 

A  -►4 

vector  r  is  expanded  on  the  special  basis  {b*'}  which  forms  the 
reciprocal  lattice. 


■*■1  *2 
=  r^b  +  r^b  ± 


(2.5) 
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Pj  =  Sj  +  kx 


(2.6) 


=  $2  +  nX  $2  =  s«d2  (2.7) 

Tj  *  y/ \  -  ir.b^  +  evaluated  in  the  (2.8) 

fourth  quadrant. 

Equations  (2.2)  -  (2.8)  have  some  features  not  found  in  the  usual 
rectangular  coordinate  plane  wave  expansion  which  we  will  now  examine. 

Our  original  lattice  provided  us  with  only  two  vectors  to  describe 
three-dimensional  space.  Therefore,  we  may  arbitrarily  choose  a  third 

-►  4.  ^  4. 

vector  d3  to  form  a  three-dimensional  basis  (di,d2,d3).  The  simplest 
choice  is  the  unit  vector 

ds  =  dj  X  d2/Bdj  x  d2»  .  (2.9) 

The  {d-j}  vectors  of  the  direct  lattice  are  also  known  as  contravariant 
vectors  and  are  a  natural  basis  for  expanding  other  contravariant 
vectors  such  as  the  positon  vector  R. 

-*■1  -*•2  ■*•3 

In  periodic  structures,  there  is  another  basis  (b  ,b  ,b  ) 
associated  with  the  direct  lattice  and  defined  by 

=  6^  (2.10) 

where  6'?  is  the  Kronecker  delta.  The  {b*^}  vectors  are  the  corresponding 
covariant  vectors  and  are  the  most  natural  basis  for  expanding 

A 

propagation  vectors  such  as  r  in  Equations  (2.5)  -  (2.8)  and  the 

A 

incident  vector  direction  s. 
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(2.11) 


) 


S  =  +  82^^  +  SjS 

(2.12) 

The  lattice  defined  by  and  b^  is  the  reciprocal  lattice  and  is 
usually  referred  to  as  k-n  space  because  of  Equations  (2.5)  -  (2,7), 

While  the  (dj,d2)  vectors  correspond  to  lattice  points,  the  (b  ,b  ) 
vectors  correspond  to  lattice  rows  as  shown  in  Figure  2.2  where  one 
reciprocal  lattice  vector  b  and  its  lattice  rows  are  given.  The  vector 
b^  has  a  direction  orthoganol  to  the  rows  and  a  magnitude  Inversely 
proportional  to  the  separation  between  the  rows.  It  is  this 
correspondence  with  lattice  rows  which  makes  the  reciprocal  basis  so 
appropriate  for  expanding  propagation  vectors. 


Figure  2.2.  Lattice  rows  corresponding  to  the  covariant  reciprocal 

lattice  vector  b  .  Each  vector  in  the  reciprocal  lattice 
has  a  corresponding  row  structure  in  the  direct  lattice. 
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Because  of  the  simple  geometry  of  the  array,  the  convention  of 
contravariant  and  covariant  vectors  is  not  really  needed.  The  important 
concept  is  the  way  the  plane  wave  spectrum  of  Equation  (2.2)  depends  on 

A 

the  reciprocal  lattice  through  r  in  Equations  (2.5)  -  (2.8). 

The  choice  of  cl^  in  Equation  (2.9)  and  the  condition  given  in 
Equation  (2.10)  allows  us  to  write 

-►1 

b  =  d2  X  d^/id^  X  d2li 

P  =  d3  X  dj/«dj  X  d2l  (2.13) 


As  an  example,  consider  the  usual  rectangular  lattice 


"l  ■  “x" 


"2  =  “z" 


s  •  SxX  +  +  s-z 


A  A 


Then  Equation  (2.9)  gives 


»di  X  d2ii  = 


The  reciprocal  lattice  becomes. 


(2.14) 


(2.15) 

(2.16) 
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A 

=  1/D^  X 

»  1/0^  z 
b^  =  -y  . 

Equations  (2,5)  -  (2.8)  then  give 


(2.17) 

(2.18) 

(2.19) 

(2.20) 

(2.21) 

(2.22) 


It  is  interesting  to  note  that  the  switch  to  a  coordinate  free 
representation  has  replaced  the  term  normally  found  in  the 
denominator  of  Equation  (2.2)  with  id^  x  d2ii  which  is  the  area  of  the 
parallelogram  formed  by  the  vectors  dj^  and  d^. 

Equation  (2.4)  for  the  vector  pattern  function  has  the  standard 
form  since  it  is  independent  of  the  lattice. 


15 


8.  THREE-DIMENSIONAL  VOLUME  ARRAYS 


Now  consider  a  volume  distribution  of  current  sources  defined  by 
the  direct  lattice  vectors  Dj,  D2  and  Dg  shown  in  Figure  2.3.  The 
reference  pont  is  chosen  to  be  at  the  origin  for  now.  The  field  from 
the  volume  array  will  be  found  by  summing  the  fields  produced  by  the 
elements  in  the  Di,D2  planes.  The  reason  for  using  the  notation  Dg 
instead  of  D3  for  the  third  lattice  vector  is  because  the  formulation 
will  be  based  on  the  expansion  of  the  previous  section  with 

di  =  Dj 


d2  =  D2 


(2.23) 


d^  =  Dj  X  Dg/JDj  X  02*  . 


*  *  *  -*•  ♦j 

The  Di,D2,  and  Dg  vectors  define  the  array  while  the  dj  and  b  vectors 
act  as  local  variables  of  the  equations.  Assume  the  currents  have  a 
plane  wave  type  of  phase  difference  but  are  otherwise  periodic. 


qmt  '  000'  qmt' 


where 


(2.24) 


(2.25) 


and  I 


(1) 

000 


(R) 


is  the  current  distribution  of  the  reference  element. 
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The  volume  array  can  be  considered  as  an  Infinity  of  parallel 
->• 

planes  with  Dj  and  D2  as  the  planar  lattice  vectors  and 


Koz  - 


(2.26) 


as  the  reference  point  of  the  tth  plane.  The  field  from  the  i=0  plane 
from  Equation  (2.2)  is 


'c’o 


00  g 


-j  6R*r^ 


2iiD,xD5i 

1  2  k=-»  n=-oo 


-*•11  C(p^^^(r^.)  x  r^.)  x  r+] 

n  .  _  _  ± 


use  for  R*(Dj  x  D2)  >  0 
-  for  R*(Dj  X  D2)  <  0  . 


The  field  from  the  plane  is. 


(2.27) 


E,(R)  =  E  fR-tDje 


-JBtD^-s 


(2.28) 


Therefore  the  total  field  is 


£=-00 


+  E^^^R) 


v  til).*  *  -JBJiSe 
♦  1  Ei‘'(R-ED  Je  ' 


(2.29) 
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-►  J 

=  S'D_  . 
e  e 


(2.30) 


We  restrict  R  to  the  region  of  validity  between  the  t=-l  and  t=l  planes; 
see  Figure  2.3.  Define 


a  =  sgn(Dg-b^) 


(2.31) 


=  sgn(Dg.(Dj  x  Dg))  . 


(2.32) 


Then  Equation  (2.29)  may  be  written  as 


-  I 


-1  Zcloe 


- j  etSe 


4=-»  2lDjxD2i 


00  00  0 

I  I  - 

k*-«>  n=-«> 


-j6(R-tDe)T(j 


[(P^^^r^)  X  r^)  X  r^] 


*  E»>(R) 


7  I 

“  Zcloe 

+  I  - 7—; — 

t=l  2lDjxD2il 


00  00  0 

I  1  - 

k=-<*>  n=-«» 


-je(R-tDe)•^o 


[(P^^^(r_  )  X  r_  )  X  r_  ] 


(2.33) 
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Rearranging  summations  gives  terms  of  the  form 


«o 


-j6t(se 

e 


(2.34) 


I 

4=1 


4-  ^ 

-j84(se  -  Oe*r.o) 
e 


(2.35) 


Using  the  properties  of  the  geometric  series 


I 


4=1 


e 


-j  0X 


(2.36) 


Equation  (2.33)  may  be  simplified  into  the  form 


E^^^R)  =  Ij^^R)  +  E^^^R)  +  E^^^R)  . 


(2.37) 


-►(1) 

Moving  the  reference  point  to  R  gives 


Er^(R)  =  — T - :r-  I  I 

2lD^  X  D2 1  |(3_oe  n 


o»  e 


-(1)  " 

-j8(R-R'  ').r^ 


-jS(r(,.De-Se) 

- [(?<'>(?,)  =<  ?„)  X 

l_g-jS(r(j»De-Se) 


(2.38) 
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>  _  _ 

2lD^  X  D2I  |(=-o<>  n=- 


m-t-  V''  /  *  r  ® 

(R)  =—:^ - 7—  I  I  - 


-jB{R-R  )T.a 


06(r_p.De-Se) 

.-- - [(?'”(?. J  X  ?.,)  X  ?.„3 


1-e 


jB(r_g«De-S0) 


(2.39) 


Ei'>(R)  « - 1 - -  I  I  - 

21D2  X  02*  k=-«o  0=-“ 


>(1)  “ 
-jB(R-R  )•% 


[(P<"(?j)  X  ?j)  X  ?j] 


(2.40) 


with 


A  +1  ■¥-^  +3 

Pjo  '  '"l'’  '"2^’  *  ®'3'> 


Equations  (2.23),  (2.13),  (2.5)  -  (2.8)  are  used  to  find  ri,r2,b. 
■^(1) 

Notice  that  Eq  (R)  does  not  depend  on  De. 


(2.41) 
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C.  IMPEDANCE  EQUATIONS  FOR  WIRE  ANTENNAS 


Place  another  antenna  element  with  transmitting  current 
{2)t  ■*•(2) 

distribution  I  (R)  and  terminal  point  R  in  the  region  of  the 
volume  array.  The  induced  voltage  in  antenna  (2)  due  to  a  unit  current 
in  antenna  array  (1)  will  be  the  mutual  impedance  given  by 


Z 


2,1 

vol ume 


(Dj,D2|Dg)  =  Z 


2,1 

+ 


{DpD2|Dg)  +  Z 


2,1 

0 


(Di,D2) 


+  z2*l(Dj,02lDg)  (2.42) 

and  is  obtained  by  combining  Equations  (1.1)  -  (1.2)  with  Equations 
(2.37)  -  (2.40).  When  I^  ^  (R)  =  I^  ^  (R)  and  R^  ^  *  R^  \  Equation 
(2.42)  yields  the  self  impedance  of  antenna  (1). 

2A  *  *  * 

Consider  first  the  term  due  to  the  plane,  Zq  (Di,D2).  If  Dj 

and  D2  are  the  initial  lattice  vectors  of  the  nonplanar  array,  then 
2 , 1 

calculation  of  Iq  (Di,D2)  leads  to  iterated  integrals  which  we  are 
trying  to  avoid.  In  order  to  find  the  nonplanar  impedance,  we  need  only 
calculate  the  impedance  once  for  the  element  in  some  two-dimensional 
lattice  in  which  it  is  totally  planar.  We  will  assume  for  this 

2  f  1  *** 

calculation  of  Iq  (0^,02)  that  and  D2  fulfills  this  condition.  In 

2,1  2,1 

Chapter  III,  we  will  use  the  Z+*  and  Z_*  terms  to  transform  the 

2,1  > 

impedance  Zq  (0^,02)  from  the  planar  lattice  to  the  original  nonplanar 
lattice. 

If  both  antenna  elements  lie  completely  within  the  plane  of  the 

2,1 

array,  then  the  self /mutual  impedance  Zq  (0^,02)  is  evaluated  using  the 
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standard  method  developed  by  Munk  [73.  In  this  procedure,  the  induced 

voltage  is  found  by  integrating  the  field  from  antenna  array  (1)  due  to 

a  filamentary  current  distribution  at  R'  over  the  filamentary 

*(2)  (Z)*3 

transmitting  current  distribution  of  antenna  (2)  at  R'  +  a  b  .  This 

point  is  one  wire  radius  a^^^  of  antenna  (2)  above  the  plane  of  the 

-*■3  .  .  . 

array  in  the  evanescent  direction  b  .  This  method  is  an  approximation 
of  Equation  (1.1)  which  has  been  shown  to  produce  excellent  results  by 
Munk  et .  al . 

A 

For  a  linear  antenna  in  the  direction  p.  Appendix  B  shows  how  the 
filamentary  pattern  factor  P^^^(r)  is  related  to  the  vector  pattern 
function  P'  '(r)  by, 

s  p  l-(r»p^^h^  )  P^^^(r)  (2.43) 


where 


P 


1 

i(i)(5(i)) 


I  i<'> 

el ement 


(t)e 


) 

dt. 


(2.44) 


J  (x)  is  the  Bessel  function  of  order  zero  and  a^^^  is  the  radius  of 

O'  ' 

A  /  1  \  A 

antenna  (1).  Equation  (2.43)  is  derived  by  averaging  pP^  ‘'(r)  around  a 

-  .  (1) 

circle  orthoganol  to  p  and  of  radius  a 

Replacing  P^^^(r)  with  p^^^P^^^(r)  in  Equation  (2.40)  and 

performing  the  integration  gives 
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-Ir 


9  1  -»■  ”^C 

z2»^Di.D2)  - 7-  I  I 


09  0 


MZ)  (2)+3  ♦(!)  * 

-j6(R  +  a  -r' 


(2.45) 


with  the  filamentary  transmitting  pattern  factor  defined  by 

^  ^(2) 

p(2)t/:x _ I -  /  T(2)t/,.  -j6t(r»p  ) 

^  ■  T(2)t,^(2).  J  ^  ®  (2.46) 

*  k''  J  element 


Equation  (2.45)  is  also  valid  if  antenna  (2)  is  outside  Region  2  of 
antenna  (1). 

2»1  *  +  2,1  +  +  + 

The  terms  Z+  (Di»D2l^e)  2.  (Di,D2|De)  are  a  straightforward 
application  of  Equation  (1.1)  and  Equations  (2.38)  -  (2.39)  and  are 
valid  for  any  element  orientation  within  the  region  of  validity 
indicated  in  Figure  2.3, 


-Zc  *  «D  e 


■*•(2)  ♦(I)  * 

-j6(R^  ^-R^ 


z!-'  =  .  i  i 

2lD,xD5l  „ 

1  2  k=-<»  n=- 


-j0(ro»De-Se) 


:r-; - t(f<»(?  )xPjxrJ.p(Z)‘(Pj 


l-e 


a'  a'  a 


(2.47) 
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r2.1  = 


2lD^xD2l 


00  00  g 

I  I 

k=-<»  n=-* 


-j0(R^  -R  ^)•^c 


where 


je(f'_a*De-Se) 


1-e 


(2.48) 


p(2)t(;:)  , 


1 


l(2)t(^(2)) 


/  I^2)^(R')  dR’  (2.49) 


is  the  vector  transmitting  pattern  function.  It  is  similarly  related  to 
the  filamentary  transmitting  pattern  factor  by 

p(2)t  ^  J(2)  (33^^^/ 1  -  (r*p^^h^)  P^^^^(r)  .  (2.50) 

Equations  (2.43)  and  (2.50)  may  be  easily  extended  to  a  piecewise 
combination  of  linear  antenna  sections. 
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D.  THE  EXCITATION  VECTOR  s 


When  using  the  volume  impedance  Equations  (2.47)  -  (2.48)  to 
calculate  the  nonplanar  self /mutual  impedance  of  an  array  defined  by 
direct  lattice  vectors  and  D2,  the  vector  is  a  dummy  vector  which 
is  usually  chosen  in  the  same  plane  as  the  antenna  elements  but  is 
otherwise  arbitrary.  Likewise,  s^  *  s*Dj  and  $2  =  s*D2  are  fixed  by  the 
scan  angle  (e.g..  Equation  (2.18))  while  s^  =  s.D^  is  arbitrary.  A 
logical  choice  for  s^  would  appear  to  be  the  one  which  makes  s  a  unit 
vector  s.  But  when  s=s,  either  (r^-s)=0  or  (r_^-s)=0  in  the  k=0,n=0 
term  in  Equations  (2.47)  -  (2.48).  This  produces  a  pole  in  the 

volume 

impedance  which  is  unacceptable  for  our  purposes.  A  safe  choice  is 
s,'0. 

The  fact  that  the  equations  blow  up  when  s  is  a  unit  vector  may  at 
first  be  disturbing  but  is  actually  related  to  the  nature  of  a  spac<» 
filled  by  a  volume  array.  The  presence  of  the  array  changes  the 
effective  constitutive  parameters  of  the  space  in  a  manner  similar  to  a 
dielectric.  The  application  of  this  principle  to  metallic  delay  lenses 
was  first  demonstrated  by  Koch  [8]. 

The  effective  vector  wave  number  k^^  of  this  artificial  medium  at  a 
given  frequency  is  determined  by  the  condition  of  a  zero  tangential 
electric  field  along  the  antenna  elements  without  an  external  voltage 
excitation.  This  corresponds  to  an  eigenfunction  of  the  periodic 
structure  and  is  equivalent  to  the  requirement 

^volume  (^0^  "  °  »  (2.51) 
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(2.52) 


With  these  conditions  fulfilled,  the  total  self  induced  voltage  is 
zero  with  a  nonzero  current  distribution  given  by  Equations  (2,24)  and 


(2.52), 


-► 

4-  4*  4-  4  4-  -jR 

qmit'  '  000'  qniJl 


4 

qm)l**^o 


(2.53) 


This  current  is  a  plane  wave  propagating  through  the  periodic  structure 

4  4 

with  phase  velocity  Uo/lk^l  and  group  velocity  doa^dlk^l.  We  would 
expect  the  presence  of  the  array  to  change  the  effective  parameters  of 

4  4 

the  medium  so  that  in  Equation  (2.52).  In  this  case,  s^^  is  not  a 

unit  vector. 

Without  the  array,  a  plane  wave  would  propagate  through  the  medium 

4  ^ 

at  frequency  with  s=s  a  unit  vector.  We  have  seen  before  how  the 
self  impedance  goes  to  infinity  when  the  volume  array  in  the  medium  is 
excited  in  this  manner.  This  gives  the  interesting  result  that  a  volume 
array  in  a  medium  with  wave  number  p  can  not  support  a  plane  wave  with 

4 

wave  number  ik^i  «  $.  The  exception  to  this  rule  occurs  when  there  is  a 

4 

null  in  the  pattern  function  in  the  direction  of  propagation  k^. 

The  periodic  volume  array  as  an  artificial  dielectric  is  not 
relevant  to  the  main  objective  of  this  report.  The  preceding  discussion 

4 

is  only  meant  to  illustrate  the  nature  of  the  vector  s  and  why  the 

4  4 

arbitrary  component  Sg=s»Dg  may  be  chosen  to  be  almost  any  value  except 

4  ^ 

the  one  which  makes  s=s  a  unit  vector.  As  mentioned  before,  $^*0  is 
usually  a  safe  choice. 
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CHAPTER  III 


THE  NONPLANAR  IMPEDANCE 

A,  GENERAL  CONSIDERATIONS 

A  planar  array  in  a  homogeneous  space  is  uniquely  defined  by  its 
♦(1) 

reference  position  R  ,  the  orientation  and  dimensions  of  the  elements 

-►(1)  -*•  ♦ 
relative  to  R  ,  and  the  lattice  vectors  Dj  and  D2.  When  the  elements 

are  tipped  out  of  the  plane  of  the  array,  the  calculation  of  the 

self/mutual  impedances  becomes  very  difficult  because  of  the  iterated 

integration.  However,  the  use  of  volume  impedances  allows  one  to 

transform  the  problem  from  the  Di,D2  lattice  to  the  03,04  lattice  where 

impedances  are  more  easily  evaluated.  In  doing  this  transformation,  the 

,  .  .  .  -►(I) 

element  s  position  R  and  orientation  remain  the  same;  only  the 
lattice  vectors  change. 

For  a  volume  distribution  of  elements  defined  by  the  direct  lattice 

4  4  4 

vectors  0|, 02,03,  the  self/mutual  impedance  has  been  derived  in  Chapter 
II  as 


.2.1 


2  1 


-  ‘0  *  zf-'tOpOalDa)  *  Z!-‘(Di,D2|03) 


.2,li 


(3.1) 
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Or  by  symmetry, 


(3.2) 

-►  ^ 

The  evanescent  direction  Dg  has  been  chosen  as  D3  in  the  first  case  and 
Di  in  the  second  case.  Combining  Equations  (3.1)  -  (3.2)  gives 


z2''(Dj.02)  .  Z^-lcOj.D,)  .  z2-'(52|Di.03) 


(3.3) 


with 


Z^’^(D2|Di-^D3)  s  zJ’^(D2.D3|Dj)  +  Z^*^(D2,D3|Dj) 


-  zJ*^DpD2lD3)  -  Zf»"(Dj.D2|D3). 


2.1 


>  >  > 


(3.4) 


Equation  (3.3)  is  the  basic  formula  for  calculating  nonplanar 
impedances.  As  long  as  the  intermediate  step  of  forming  a  volume  array 
is  valid,  the  impedance  terms  in  Equation  (3.4)  are  always  calculable 
and  have  a  strong  exponential  convergence. 

If  the  array  is  still  nonplanar  in  the  02,03  lattice  but  planar  in 
the  03,04  lattice,  another  transform  may  be  possible. 


(02,03)  =  Z 


2,1 

0 


(03,04)  +  Z 


2.1 

A 


(03102-^04)  . 


(3.5) 


Combining  with  Equation  (3.3)  gives 
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The  same  vector  pattern  functions  P^^^r),  and  filamentary 

pattern  factors  P'  '(r),  P'  (r)  are  used  throughout  these  transforms 
since  the  elements  maintain  the  same  orientation. 


B.  DIPOLE  ARRAY 


Consider  an  array  of  linear  dipole  antennas  with  arbitrary 

A 

orientation  p  in  a  rectangular  grid  with  Interelement  spacings  and 
as  shown  in  Figures  1.1  and  3.1.  Equation  (3.3)  will  be  used  to  find 
the  self  impedance.  For  this  geometry 


Figure  3.1.  Single  dipole  used  in  the  example. 


We  will  assume  a  linear  phase  difference  between  elements  characterized 
■*.  -jes'R 

by  a  plane  wave  Ege  excitation; 


(3.8) 


and  therefore, 


s 

s 


(3.9) 


The  antenna  elements  have  a  half  length  i  and  radius  a.  We  will 
choose. 


Dg  =  2.5  £  p 

S3  =  0  .  (3.10) 

This  choice  Insures  a  valid  volume  array  and  insures  the  dipole  array 
will  be  totally  planar  in  the  ^2*^3  lattice  (z-p  plane).  A  larger 
means  a  faster  converging  2^(0^,02103)  term  but  a  slower  converging 
^o^^2*^3^  term.  Equation  (2.45)  is  now  used  to  calculate  the  planar 
20(02,83)  impedance  and  Equations  (2.47),  (2,48)  and  (3.4)  are  used  to 
find  the  difference  2^(02|0j-^03)  impedance.  Adding  gives  the  self 
impedance  of  the  nonplanar  dipole  array.  Figure  3.3  shows  the  scan 
impedance  for  the  case. 
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PLANE  OF  INCIDENCE 


Figure  3.2.  Incident  angle  parameters  a  and  n. 


Dx  =  =  1  cm 

A 

P  =  -y 
i  =  .68  cm 
a  =  .05  cm 
e  =  2.17  . 

The  incident  angle  parameters  a  and  n  are  illustrated  in  Figure  3.2. 

As  another  simple  example,  the  mutual  impedance  is  calculated  for 
two  y  directed  dipole  arrays  having  the  same  parameters  as  the  first 
example  as  a  function  of  their  separation  in  the  z-y  plane.  The 
positions  of  the  two  reference  elements  are 


IMflG 

-300  -200  -100  „  0  100  200  300 


SCRN  IMPEDANCE 


BLPHfi  =0.00 
F  (GHZ)  =11.  00  o 

F  (GHZ)  =5.  00  A 

F  (GHZ) =6. 00  + 

F  (GHZ)  =7.  00  X 

F  (GHZ) =8. 00  A 


Figure  3.3.  Scan  impedance  for  an  array  of  y  directed  dipoles  for 
various  angles  of  incidence  and  frequencies. 
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=  origin 


r(1) 

*  Oyy  +  O^z 


(3.11) 


where  the  offsets  Oy  and  O2  are  measured  in  cm.  Figure  3.4  shows  the 
geometry  of  the  problem  and  figures  3.5  through  3.9  show  the  scan 
impedance  for  various  offsets  and  angles  of  incidence  with  a*90°.  In 
this  case,  D3  was  chosen  as 


which  was  sufficient  for  the  examples  treated. 


C.  LOOP  ARRAY 

A  common  method  for  constructing  nonplanar  phased  arrays  is  to  etch 

one  row  of  elements  on  printed  circuit  boards  as  in  Figure  3.10  and  then 

stack  the  boards  to  form  the  Di,D2  planar  lattice  shown  in  Figure  3.11. 

By  choosing  D3  in  the  plane  of  the  pc  boards,  the  0^,03  plane  may  be 

2,1  -*• 

used  to  calculate  Zg  (01,03)  for  all  self  and  mutual  impedances  between 

antenna  elements.  As  before,  this  planar  impedance  is  readily 

2,1  -►  + 

transformed  into  Zg  (01,03)  for  the  nonplanar  lattice.  This  is  perhaps 
the  most  useful  application  of  this  method  because  of  the  ease  with 
which  complex  current  distributions  and  modes  may  be  intermixed. 
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fiLPHR  =90.00  0Z=  0.20 

p(GHZ)=il.0C  o  0T=  0.00 

F (GHZ) =5. 00  A 
F  (GHZ) =6. 00  + 

F (GHZ) =7. 00  X 
F  (GHZ)  =8.  00  c> 


Figure  3.5.  Mutual  impedance  between  y  directed  dipole  arrays  with 
separation  0^=0. 0,  0^=0. 2. 
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SCAN  IMPEDANCE 


RLPHA  =90.00 
FtGHZ)=ii.OO  o 
FCGHZ)=5.00  A 
F(GHZ:=6.00  + 

F  (GHZ) =7.00  X 

F  (GHZ) =8.00  ^ 


Figure  3.6.  Mutual  impedance  between  y  directed  dipole  arrays  with 
separation  0^*0. 0,  0^=0. 35. 


0Z=  0.35 
0T=  0.00 
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-300  -200  -100 


SCfiN  IMPEDRNCE 


o 


CO 

I 


RLPHB  =90.00  0Z=  0.50 

F  (GHZ) =4. 00  o  0T=  -0.50 

FtGHZl=5.00  A 
FCGHZ1=5.00  + 

F  (GHZ) =7.00  X 
F  (GHZ)  =6.  00 


Figure  3.8.  Mutual  impedance  between  y  directed  dipole  arrays  with 
separation  0y»-0.5,  0^*0. 5. 
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TERMINALS  FOR  PC  BOARD 

FEED  NETWORK 


Figure  3.10.  Single  row  of  loop  antennas  etched  on  a  copper  clad  pc 

board  with  interelement  spacing  Rotation  angle  6  is 

used  in  the  example. 


Figure  3.11.  Nonplanar  array  formed  by  stacking  pc  boards. 


We  will  consider  as  an  example  an  array  of  square  loop  antennas 
arranged  in  a  hexagonal  lattice  in  the  x-z  plane.  The  loops  will  be 
rotated  in  an  angle  6  in  the  plane  of  the  pc  boards  as  shown  in  Figure 
3,10,  The  pc  boards  will  be  stacked  so  that  they  are  rotated  by  an 
angle  <|>  relative  to  the  plane  of  the  array.  Figure  3.12  is  an  edge  on 
view  of  the  pc  boards  and  the  plane  of  the  array  with  Dj  pointing  into 
the  page.  The  hexagonal  lattice  with  e=0  and  <^=0  is  shown  in  Figure 
3,13;  only  the  edge  of  the  loops  are  visible.  By  definition,  «Dji  * 

>D2i  and  the  angel  between  Dj  and  D2  is  Tr/3  for  this  lattice.  An 
examination  of  the  reciprocal  lattice  would  show  that  it  is  also 
hexagonal , 

We  will  assume  the  same  symmetrical  cosine  current  distribution  for 
transmitting  and  receiving  as  used  by  Kent. 

The  planar  self  imipedance  will  be  evaluated  in  the  (Dj.Dj)  plane 
and  then  transformed  to  the  (Dj.Dg)  plane.  A  convenient  choice  for  D3 
is 


Dj  =  2)i(sin(})X  -  costy) 


(3.13) 


and  s^  is  chosen  to  be  zero  as  usual. 

Figures  3.14  -  3.17  show  the  self  impedance  of  this  array  in  a 
homogeneous  e=2.17  space  with  the  following  parameters: 


iDin  =  1,07  cm 
t  =  .68  cm 
a  =  ,05  cm 
E  =  2.17  cm 
e  =  0“,  45" 

^  =  0°,  45°  . 
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Figure  3.12.  Stacked  pc  boards  seen  edge  on  and  rotated  by  the  angle 

♦  ♦ 


A 

z 


Figure  3.13.  With  9=0  and  (|»=0,  only  the  edge  of  the  loops  can  be  seen 
looking  down  onto  the  hexagonal  lattice. 
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o  * 

o 

' 

RLPHR  =0.00 

THETPI=  0.00 

FtGHZ)  =4.  00 

O 

PHI=  0.00 

FrOHZl  =5.00 

FCGHZJ  =6.  00 

+ 

FCGHZ)  =7.00 

X 

FCGHZ)  =8.00 

o 

Figure  3.14.  Self  impedance  for  a  loop  array  with  9=0°  and  ((»=0°. 


O  J 
o  -J 
o 

I 

RLPHR  =0.00 
FCGHZJ=1A.00  o 
F[GH2J=5,00  A 
FtGHZ)=6.00  + 

F  [GHZ)  =7.00  X 
FCGHZ)=8.00  o 


THETfl=  0.00 
PHI=  45.00 


Figure  3. 16,  Self  impedance  for  a  loop  array  with  0=0®  and  4>=45®. 


CHAPTER  IV 


CONCLUSION 


The  method  of  transforming  the  Impedance  between  different  lattices 
as  presented  in  the  previous  chapters  can  be  a  very  useful  tool  for 
analyzing  nonplanar  arrays.  Its  value  depends  heavily  on  the 
application.  For  unusual  current  distributions  with  several  overlapping 
modes  all  confined  to  the  same  plane  outside  the  plane  of  the  array, 
this  method  would  be  the  only  practical  approach.  However,  for  elements 
with  a  three-dimensional  structure,  the  iterated  integral  approach  may 
provide  the  only  solution.  Either  way,  the  researcher  must  pay  careful 
attention  to  the  geometry  of  the  problem.  Given  a  choice,  the  transform 
method  can  provide  an  exponentially  convergng  double  summation  in  all 
its  terms  as  opposed  to  the  algebraic  convergence  obtained  using 
iterated  integrals. 

Only  analysis  has  been  considered  so  far.  Equally  important  for  a 
model  is  its  design  potential  or  the  ease  with  which  the  equations  can 
be  broken  down  into  simpler,  more  manageable  form.  Both  the  transform 
method  because  of  its  dummy  lattice  vectors  and  the  interated  integral 
method  because  of  its  cumbersome  complexity  fail  miserably.  The  power 
of  the  plane  wave  expansion  comes  from  the  simple  cause  and  effect  form 
which  it  takes  for  a  totally  planar  array.  Future  research  in  nonplanar 
arrays  should  address  this  issue. 
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APPENDIX  A 


DERIVATION  OF  THE  FIELDS  FROM  A  PLANAR  ARRAY 


In  this  appendix,  the  field  from  an  infinite  planar  lattice  of 

current  sources  will  be  derived.  The  geometry  is  shown  in  Figure  2.1 
•►(1) 

with  R  =  origin.  The  lattice  is  defined  by  the  two  vectors  d^  and  d2 
and  the  current  Iqni(R)  associated  with  the  q,m^"  lattice  point  Rq^  will 
be 


00  qm 


(Al) 


with 

R„„  *  (^2) 

qm  ^  1  z 

-j  0s  •R 

upon  excitation  by  a  plane  wave  E^e  .  Equation  (Al)  is  simply  a 

statement  of  Floquet's  theorem.  Alternatively,  this  current  may  be 
impressed  on  the  elements  by  an  external  generator. 
current  distribution  of  the  q=0,  m=0  element  relative  to  the  origin. 
The  total  current  is 


I{R)  = 


00 

I 


qs-00 


I 


I,™*"' 


1 1 


00 


(R- 


%m)® 


-j  8S  ‘R 


qm 


qs-oo 


(A3) 


(A4) 
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I(R) 


4-  -► 


■  I  W"')  I  i 


-  j  Ss  ‘RqiT, 


element 


6{R-Rqn,-R')}dR' 


q=-»  m=-co 


(A5) 


using  the  properties  of  the  delta  function. 

The  scattered  field  will  be  found  by  solving  the  differential 
equation  for  the  vector  potential  A(R): 

V^A  +  rA  =  -  yI(R) 


2  2 
0  *  (I)  we  , 


The  problem  is  simplified  by  considering  the  scalar  equation. 


(v^+0^)G(r)  =  w{  I  I  e''^^^’'^^"'«(R"-Rqn,)}  • 
q=-“  m=-'» 


Integrating  both  sides  of  this  equation  by  /Iqq(R')  with 


R"  =  R-R' 


(A6) 


(A7) 


(A8) 


gives 


A(R)  =  J  lQp(R')G(R-R')dR' 
element 


where  6(R")  is  the  scalar  Green's  function. 


(A9) 


(AlO) 
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Evaluating  G(R) 


The  primes  in  Equation  (A8)  will  be  dropped.  Let  the  Fourier 
transform  of  G(R)  be  defined  by 

F  (I)  =  ///  e"^®^’S(R)dR  .  (All) 


Applying  the  transform  to  Equation  (A8)  and  substituting  Equation  (A2) 
gives 


S^(i4-E)F(5)  = 

-  -  -jS(qsi+mS2)  jB(q?i+m52)i 

-ui  2,  I  e  e  } 

q=_o*  ni=-«» 

(A12) 

where 

A 

5j  •  S.dj 

Sj  •  S.dj 

(A13) 

El  -  E-di 

Using  the  identify 

^2  =  C'd2  . 

-  Jnui„t  2.  -  2. 

I  e  =:r  ^  ) 


n=-oo 


n=-«B 


(A14) 


with  the  substitutions 


I  =  X  I  6{Ci-{Sj+kX)) 


q=. 


?  jSm{C2-S2)  V  / 

1  e  *  X  ^  6(52-(s2+nX)) 


in=-oo 


n  =  <-eo 


(A16) 


Substituting  Equations  (A16)  into  Equation  (A12)  we  obtain 


,  .„x^  -  .  «(Ei-(Si-H<»  «(E2-(S2+nx)) 

f(E)  •—  I  I  - 


k=-»  n=-« 


(A17) 


Applying  the  inverse  transform  to  Equation  (A17): 


-Xvi  »  00 

G(l*)  ii/  I 

-00  k®-* 


T  -jBR*5 

I  e 

n=-oo 


6{Cj-(Sj+kX)6(C2-(S2+nX)) 

l-l^K 


-► 

dK 


(A18) 

So  far,  only  two  vector  directions  have  been  chosen  dj  and  d2.  A 
third  independent  vector  is  needed  to  completely  specify  any  direction 
in  space  as  a  linear  combination  of  these  three  vectors.  We  will 
choose 

da  =  dj  X  d2  /  idj  X  d2l  .  (A19) 

This  choice  provides  the  greatest  simplicity  in  the  resulting  equations. 

> 

With  the  three  (contravariant)  vectors  dj,d2,da  specified,  we  can  form 

-►I  -►O  -fO 

the  corresponding  covectors  b  ,b  ,b  according  to  the  formula 
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=  6^’  (A20) 

6^  =  Kronecker  delta. 

Equations  (A19)  and  (A20)  yield  the  following  equivalent  forms: 

♦  I  4-  -►  4-  -► 

=  dg  X  dj  /  idj  X  ^2* 

■*■0  ■¥•*■■*■  * 

*  d2  X  dj  /  »dj  X  d2* 

P  *  dg  .  (A21) 

Now  form  a  transformation  of  the  integral  to  the  system  of 

coordinates  where 


*  *9  *2 

K  =  C^b'  +  +  Cjb 


^i  '  ^-^i  • 


(A22) 

(A23) 


This  definition  of  Ci  is  consistent  with  Equations  (A13).  Under  this 
transform.  Equation  (A18)  becomes. 


-Xu 


00  00 


G(R) 


I  I  JJJ  e 


-jS(RUl+R2S2+f^^^3) 


4.  lajxJj'  k...  .. 


5{5i-(Sj+kx))6{52-S2+nX)) 


1  - 


(A24) 
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G{R)-— ^  I  I  ^-JB(Rlri.R2r2) 

AirZidjXd^ii  k=-«  n=-« 


-  (l-.r,b'*r/.2)  “'3 


(A25) 


with 


-*•  1+  04-  -3^ 

R  =  R^dj  +  R'^d^  +  R  dj 


i  ■*■  '*‘i 
R  =R*b^ 


(A26) 


Pj  =  Sj  +  kA 


r2  =  $2  +  nX  . 


(A27) 


The  integral  in  Equation  (A25)  is  a  standard  integral  often  found  in 
oscillating  mechanical  systems  and  can  be  easily  evaluated  by  residues. 


3 

-jBR  53 


i  2 

53  -  (l-l5jb^+52b^i^) 


-irj  ^jBR^r^ 
tl..  Z2.2^  ® 


(A28) 


3  - 

use  -  if  R  =  R»b  >  0 
3  * 

+  if  R  =  R*b  <  0 


and 


/ *Z  2 

Pa  =  V  1  -  iPib  +r2b  II  evaluated  in  4^*^  quadrant.  (A29) 
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Substituting  Equation  (A28)  into  Equation  (A25)  we  obtain. 


G(R)  = 


-jXu 

4Trl2[jx52l 


-jeR*r+ 

00  00  e 


(A30) 


r 


+ 


+ 


r  it3 


(A31) 


Equations  (AlO)  and  (A30)  now  combine  to  give  us  the  vector  potential 
A(R)  for  the  planar  array: 


A(R)  = 


-jXuI(o) 
4ir  13^x5211 


00  00 

I  I 

n=-»  k=-» 


-j BR«r+ 
e 


^3 


•¥  A 

P(r^) 


(A32) 


P(r) 


1 


-►  j$r»R'  -► 

i  I(R')e  dR'  . 

elements 


(A33) 


The  subscripts  have  been  dropped  from  the  current  and  I(o)  refers  to  the 
magnitude  of  the  reference  current  at  the  origin;  this  is  usually  the 
terminal  current.  Equation  (A33)  is  the  definition  of  the  vector 
pattern  function  P(r). 


The  Electric  and  Magnetic  Fields 

The  magnetic  field  is  derived  from  Equation  (A32)  by 

H  =  V  X  A  (A34) 
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with  the  vector  identity 


V  X  (pf)  *  vf  X  p  (p  is  a  constant  vector) 


(A35) 


which  gives 


H(R) 


1(0) 


m  ee  0 

I  I 


-jSr+«R 


2ld,xd„il  . 

1  2  0=-*  k=- 


(P(rJ  X  rj  .  (A36) 


Applying  Ampere's  equation 


E  =  —  7  X  H 
Jue 


(A37) 


and  Equation  (A35)  to  Equation  (A36)  gives 


A 

-j6r+*R 


•♦■A  ^  A 


+  >  Zcl(o)  *  »  e 

E(R)  =  — r~r“  I  I  ~  [(P(r+)  X  r^.)  x  r^.]  (A38) 

2Bd,xd,ii  ^  3  -  -  - 

1  2  n=-®  k--* 


with. 


Zr  =  /u/e  . 


{A39) 


-(1) 

when  the  reference  point  is  located  at  R  , 


ZrI(R^^^)  *  we 

E(R)  =-T-; —  I  I  - 

2id,xd«B  „  . 

1  2  n=-«>  k=-'» 


"  -  ■*•(!) 
•j6r+*{R-R  ) 


[(P{r,)  X  ;,)  X 

(A40) 


P^^^(r)  =-T7Tr/  I(R')e^  dP'; 


I(R^^)  element 


(A41) 


■^(1) 

R'  is  a  vector  from  R  to  the  integration  point. 
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APPENDIX  B 


PATTERN  FACTORS  FOR  DIPOLES  AND  LOOPS 


In  this  appendix,  pattern  factors,  vector  pattern  functions,  and 
current  distributions  used  in  the  examples  will  be  considered. 

For  the  case  of  a  linear  round  wire  with  a  current  distribution 
I(jl')  and  direction  p  shown  in  Figure  B.l,  the  vector  pattern  function 
is  defined  to  be, 

A  ♦  * 

-►aP  j6R**r> 

P(r)  =  dr  (Bl) 

element 


from  Equation  (2.4).  Therefore, 


P(r) 


I^(o) 


1 

27  ^ 

o 


2ir  f 

J  ^ 


/ 

j$R  -r 


dt'  d((i' 


(B2) 


R'  =  i'p  +  a(tj  coS())  +  tg  sin(|»)  (B3) 

where  {p,tj,t2)  form  an  orthonormal  basis,  a  is  the  wire  radius,  and  I^ 
stands  for  the  filamentary  current  distribution.  Separating  the 
integrals  gives 

^  .  1  2ir  jea{r.  cos*  +  r.  sin*) 

P(r)  »  pP(r)  I  e  ^2  d*}  (B4) 

0 
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Figure  B.l.  A  linear  round  wire  with  a  current  distribution  pl()i'). 


r  =  r.t2 


and 


Tp  =  r-p 


(BB) 


P(r)  «  -7^  I^(l')e'^®^''*Pdt- 

1^0)  , 

*1 


(B6) 


is  the  filamentary  pattern  factor. 

It  can  be  shown  from  the  properties  of  the  Bessel  function  that 


O 


{B7) 
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Therefore, 


P{r)  =  ?P(P)  hi  ^ 


Since  r  is  a  unit  vector. 


r,.2  +  r^.2  +  p  2  s  1 

ti  t2  p 


and  so 

P{r)  *  pP{r)jQ($a  /l-rp2) 
which  is  Equation  (2.43).  Similarly, 

P(r)  =  pp^(r)jQ($a  yrr^) 


(B8) 


(B9) 


(BIO) 


(BID 


Dipole  Current 

The  filamentary  current  distribution  on  a  dipole  under  scattering 
conditions  is  given  by 

I^(t')  =  (A  +  Bcos(Bt'))  U'|<4  (B12) 

illustrated  in  Figure  B.2.  The  constants  A  and  B  are  determined  by  the 
conditions 

1(0)  =  1 

■  0 

=  t  +  At  (B13) 

e 
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Figure  B.2.  Scattering  (cosine)  current  distribution  on  a  dipole. 


in  which  At  is  a  correction  due  to  end  effect  capacitances  and  is  set 
equal  to  zero  in  the  examples  of  Chapter  III.  Equations  (B12)  and  {B13) 
yield, 

cos(Btg) 

^  ~  cos($ig)-l 

1 

6  ■  r-lSsTeV)  •  (B‘5) 

e 

The  scattering  current  pattern  factor  is  then  found  to  be, 
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} 


P(rp)  =  2Ai 


sin(erpt) 

V 


2B 


rpCos(BJt)s1n(erpt) 


sin(g£)cos(erpt) 


(B16) 


This  expression  combined  with  Equation  (BIO)  gives  the  vector  pattern 
function  for  a  scattering  current  distribution. 

The  filamentary  current  distribution  on  a  dipole  under  transmitting 
conditions  is, 


I^(£')  =  sin(B(£g-|^'l))  for  |£'|<£  (B17) 

illustrated  in  Figur*'  B.3.  The  transmitting  pattern  factor  is  obtained 
using  Equation  (2.46)  and  yields. 


P^{r„)  = - ; -  {-cos(B£-)  +  cos(6A£)cos(Br^£) 


-  rpSin( BA£)sin(Brp£)}  . 


(B18) 


The  vector  transmitting  pattern  function  is  found  from  Equation  (Bll). 


61 


0 


Figure  B.3.  Transmitting  (sine)  current  distribution  on  a  dipole. 


Loop  Current 


The  cosine  current  distribution  used  by  Kent  for  a  rectangular  loop 
is 


^  cos (at') 

'  *ros28(>*b)  (819) 

and  is  shown  in  Figure  B.4.  Each  side  contributes  the  following  terms 
to  the  pattern  factor: 

-2 

^os(ac)(l-rp2)  {s1n(Bc)-cos(Br^jb)sin8(c-b) 

-  r^sin(Brjjb)cos(B(c-b))}  (B20) 
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Figure  B.4.  Rectangular  loop  antenna. 


.  2e'Ba2ra 

Scos($c){l-r^j2)  {sin(Bb)cos(Brjjb)  -  rjjSin(Brjjb)cos(Bb)} 

(B21) 

*  ..  2sin(Bbrb) 

P2(rHP4(r)  •  Bcos(6c)(l-r  2)  ^-'•a“s(6(c-b))  -  Js1n(8(c-b)) 

d 

+  r  e'^*^^^^®cos{Bb)  +  j  e'^^^^''®sin{6b)}  (B22) 

a 


where 


The  vector  pattern  function  is  therefore. 


P(?)  •  {Pj  o„(ea/r:?^)[P2(r)  +  P,(p)] 

t  0„(Sa/Tn^[P,(r)  ♦  P3(0]t  . 
The  transmitting  vector  pattern  function  is 

“►4-  A  -►  A 

P^(r)  =  P(-r)  . 


(B24) 


(B25) 
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